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L  INTROOrCTIOX 


During  the  period  from  July  1,  19SG  ihrougli  August  31.  lOjD.  tlie  execution  of  contract 
No.  F19G_'>-.^G-C-0Ui9  was  diiccted  toward  the  development  and  implementation  of  numerical 
ntethods  for  simulating  and  describing  the  transport  of  electrons  and  photons  in  lav  ered  une- 
dimensional  and  multiregion  two-  and  three-dimensional  structures  such  as  metallized 
semiconductors.  In  addition  to  microelectronic  device  structures,  the  numerical  methods 
d.;vo!oped  were  aDo  applied  to  liquid  media  (i.c.  H^Oj.  The  tran.sport  simulations  were  used 
to  calculate  charge  distributions  and  to  predict  dose  distributions  in  these  media,  resulting 
from  electron  beams  and  sources  of  x-  and  gamma  irradiation.  .A.nothei  important  aspect  of 
the  electron  transport  calculation  activity  during  this  contract  was  the  extension  of  the  range 
of  our  simulation  capability  to  electron-phonon  scattering.  The  effects  of  this  interaction  and 
the  effect  of  the  presence  of  applied  electric  fieids  on  the  energy  spectra  and  spatial  distribution 
of  electrons  in  SiOs  were  examined,  with  particular  attention  to  the  promotion  of  electron 
energies  from  the  subionization  region  to  energies  above  the  ionization  threshold  and  the 
resulting  electron  multiplication.  Some  of  the  results  of  the  electron  and  photon  transport 
calculations  are  reported  in  seven  technical  papers  in  which  the  principal  investigator  shared 
authorsnip. 

The  electron  and  photon  transport  calculations  were  made  using  the  method  of  dis  rete 

ordinates  and  the  .Monte  Carlo  trajectory  simulation  method.  The  discrete  ordinates 

calculations  were  made  using  the  ONETRAN^*^  computer  code.  With  this  code  we  calculated: 

(1)  charge  distributions  and  energy  deposition  in  nine  materials  due  to  the  incidence  of 

monodireciionai  electron  sources  for  several  energies  and  incident  angles;  (2)  electron  and 

ohoton  transport  in  H^O;  (3)  energy  deposition  in  and  electron  transmission  energy  spectra 

from  A1  for  the  purpose  of  comparison  with  independent  transport  calculations  made  using  the 

method  of  streaming  rays.  The  Monte  Carlo  method  was  used  to  perform  electron  transport 

calculations  in  one-,  two-  and  three-dimensions.  Most  of  the  Monte  Carlo  calculations 

[21 

involved  the  use  of  the  ITS  code  series.  However,  some  of  the  one-dimensional  calculations 
were  performed  with  Monte  Carlo  programs  which  we  wrote  expressly  for  the  purpose:  in 
particular  an  electron  transport  calculation  to  test  a  new  scattering  cross  section  formulation: 
and  an  electron-phonon  scattering  simulation. 
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This  report  is  organized  into  five  sections  in  addition  to  this,  the  introduction.  Section  II 
is  a  discussion  and  presentation  of  results  of  one-diniensional  discrete  ordinates  calculations  of 
electron  and  photon  transport.  Section  III  presents  the  results  of  electron  transport  Monte 
Carlo  calculations  in  one-,  two-  and  three-dimensions  and  shows  comparisons  of  these  results 
with  transport  solutions  obtained  by  other  methods.  Section  IV  consists  of  a  discussion  of 
electron-phonon  scattering  and  some  of  the  results  obtained  with  our  Monte  Carlo  simulation 
of  electron  transport  in  the  1  to  20  eV  energy  range.  Section  V  is  a  discussion  of  the 
formulation  of  and  progress  made  towards  a  discrete  ordinates  solution  to  the  electron 
transport  problem  in  the  1  to  20  eV  energj-  range.  The  final  section,  VI,  is  a  list  of  the 
references  cited  throughout  this  report. 


IL  DISCRETE  ORDINATES  C AI  CCLaTIONS  .Q£  ELECTRON  AND  PHOTON- 
TRANSPORT  1^  ONE-DIMENSION 
1 .  Introdiiction 


Prior  to  tliL-  start  of  tliis  reporting  period,  we  had  acquired  a  substantial  aiiiounl  of 
experience  iii  adapting  the  multigroup,  discrete  ordinates  method  to  the  calculation  of  electron 
energy  deposition  profiles  a.id  electron  flux  distributions  in  several  materials.  In  view  of  this, 
it  was  deemed  appropriate,  during  the  performance  of  this  contract,  to  perform  a  systematic 
study  of  the  range  of  validity  of  this  method.  The  conclusions  reached  as  a  result  of 
performing  this  work  were:  1)  that  the  simple  continuous-slowing-down  approximation 
(CSD.\)  electron  transport  model  was  adequate  for  electron  transport  calculations  in  the 
energy  range  100  keV  to  ~1.5  Me\'  for  a  wide  range  of  scattering  materials;  and  2)  above  this 
energy  range,  the  CSD.A  pro\ed  inadequate.  This  second  conclusion  led  to  an  extensive 
amount  of  effort  to  account  for  radiative  energy  losses  and  knock-on  secondary  electron 
production  in  our  transport  calculations.  To  account  for  these  two  electron  energy  loss 
mechanisms,  we  began  an  effort  to  apply  the  discrete  ordinates  method  to  electron-photon 
transport.  We  shall  first  discuss  the  results  of  the  application  of  the  method  of  discrete 
ordinates,  incorporating  only  the  CSDA  model,  to  the  determination  of  electron  energy 
deposition  and  charge  distributions  in  several  materials.  Following  this  will  be  a  discussion  of 
our  application  of  the  discrete  ordinates  mathod  to  electron-photon  transport. 

Anolication  of  the  Discrete  Ordinates  Method  to  Electron  Transport  hi  the  0,1  •  VA  MeV 
Energy  Ranee 

The  multigroup  discrete  ordinates  code  ONETRAN^*^  was  used  to  calculate  energy 

deposition  profiles  and  charge  distributions  in  a  number  of  solids  due  to  the  incidence  of 

electron  beams  ranging  in  energy  from  0.1  to  1.4  MeV.  In  these  calculations  electron  energy 

loss  was  allowed  only  through  application  of  the  continuous-slowing-down  approximation 

(CSDA).  Thus  energy  downscatter  was  allowed  only  between  adjacent  energy  groups.  Energy 

losses  through  radiative  processes  or  through  catastrophic  (knock-on)  collisions  were  not 

considered.  The  primary  motivation  for  this  study  was  to  determine  the  range  of  validity  of 

the  method  of  discrete  ordinates  under  the  above  straightforward  assumptions.  The  transport 

calculations  were  made  for  several  materials  over  a  wide  range  of  atomic  numbers  (Z=4  -  92). 

A  large  body  of  data  '  ,  experimental  measurements  and  Monte  Carlo  calculations  of  energy 

deposition,  was  available  to  us  for  comparison  purposes.  Some  of  the  results  obtained  will  be 
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presented  here.  However,  the  complete  body  of  results  and  coinpaiisons.  consisting  of  3  flux 
plots.  14  cnerg>  deposition  plots,  a  transmission  spectrum  and  a  reflection  spectrum,  is 
presented  in  Ref.  5. 

Since  the  basic  principles  of  i!ie  discrete  ordinates  method  are  presojite-l  brielly  in  Ref.  5 
and  extensively  ei.sewhere  le.g.  Ref.G).  it  is  sufficient  to  state  here  that  in  the  inultigroup 
approximation,  the  electron  energy  range  is  divided  into  a  number  of  energy  groups  of  width 
AE^:  the  nionoenergeiic  traiisport  eciuation  is  then  solved  for  each  encrgt  group  flux.  Oj-  with 
the  removal  cro.s.s  section  Gaj  for  group  g  given  by  CSDA  as 

(Tag  —  S(Ej)/Ar..j,  (1) 

wliere  S(Es)  is  the  collision  stopping  power  for  electrons  with  energy  Eg.  If  the  group 
structure  is  arranged  s.ich  that  increasing  group  index  g  correspond.s  to  decreasing  energy,  then 
the  electrons  removed  from  group  g,  OagOg.  become  the  inscatter  source  to  group  g+1.  The 
group  flux,  0g,  i.s  a  function  of  po.sition  and  angle.  The  angular  dependence  is  treated  by 
resolving  the  flux  into  a  number  of  angular  components,  each  corresponding  to  a  particular 
discrete  ordinate,  and  then  solving  the  transport  equation  for  each  of  these  components.  The 
spatial(x)  dependence  of  the  flux  is  solved  for  by  the  use  of  the  diamond  difference 
approximation,  or  in  the  case  of  O.NETRAN,  the  linear  discontinuous  finite  element 
approximation. 


If  we  define  the  angular  flux  along  a  particular  discrete  ordinate,  say  p,,  as  ^,j(x).  then 
we  can,  by  integration  over  angle,  obtain  the  group  scalar  flux  C>g{x).  if  the  integration  is 
performed  numerically,  tiien 

M 


(2) 


i  =  l 


where  the  w,  are  the  numerical  quadrature  weights  corresponding  to  the  discrete  ordinate 
direction  p,.  A  sample  plot  of  ^gix)  is  shown  in  Fig.  1  for  the  case  of  a  1  MeV 
monodirectional  electron  source,  incident  direction  Hq  =  1.0,  located  at  x  =  0  in  an  infinite 
medium  of  Be(Z=:4).  The  ONETRAN  calculation  was  made  for  200  energy  groups  of  equal 
width.  Shown  in  Fig.  1  are  the  group  scalar  fluxes,  as  a  function  of  distance,  for  every  tenth 
group,  beginning  with  a  top  energy  of  0.95  MeV.  From  a  figure  such  as  this  one  can  gain  an 
appreciation  for  the  change  in  the  spatial  characteristics  of  the  charge  distribution  as  the 
electron  energy  decreases.  In  Ref.  5  other  scalar  flux  plots  are  given  for  1  MeV  electron  in 
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higher  Z  materials.  Fe'Z  =  26)  and  Ta(Z=73). 


Depth  (g/cfn2) 


Figure  1.  Electron  group  scalar  flux,  3.3  calculated  with  ONETRAN^*^  in  an  infinite  medium 
of  Be  for  1  MeV  electrons  incident  at  an  angle  of  O'  (Ref.  5)  . 

From  these  it  can  be  observed  that  the  directionality  of  the  source  electrons  is  sooner 
‘'forgotten”  in  the  high  Z  materials  than  in  materials  of  lower  Z.  Since  the  scattering  in  low  Z 
materials  is  more  anisotropic  than  for  high  Z,  the  electrons  approach  diffusive  behavior  sooner 
(at  higher  energy)  than  in  low  Z  materials. 
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The  energy  deposition  profile  Wix;  is  calculated  as 


G 

\V(xj  =  Yi  ’?v(x)  S(Ey)  +  c)o(x)  S(EJ  .  (3) 

where  Oo(x)  is  the  unscattered  flux  at  the  source  energy  E^.  and  G  is  the  total  nurnber  of 
groups.  Typical  energy  deposition  profile  plots  are  shown  in  Figs.  2,  3  and  -4.  In  all  three  of 
these  figures  we  compare  our  O.NETRAN  result  (solid  curve;  with  the  .Monte  Carlo  calculation 
(histogram)  and  c.xperimcntal  calorimetric  measurements(circles)  of  Lockwood,  Miller  and 
Halbleib*  .  In  Fig.  2,  the  energy  deposition  is  plotted  for  a  0.1  .MeV  electron  beam  norm  ally- 
incident  on  Be.  The  agreement  with  Monte  Carlo  is  good,  but  not  perfect.  The  agreement 
with  experiment  is  not  very  good,  but  then  neither  is  the  agreement  between  experiment  and 
Monte  Carlo.  Fig.  3  shows  the  energy  deposition  results  for  a  slant  beam  (60*)  of  0..5  Me\' 
electrons  incident  on  .Mo.  The  overall  agreement  among  the  three  metiiods  is  much  bettm-  in 
this  case.  The  results  of  a  OiNET'IAN  calculation  for  a  three-layered  scattering  medium 
(.■\1/Au/A!)  is  shown  in  Fig.  4.  A  1.0  Me\'  electron  beam  is  normally  incident  on  the  vacuum 
boundary  of  the  leftmost  layer  (.-M).  As  can  be  seen,  the  agreement  is  very  good. 

2.  Extension  of  the  Discrete  Ordinate.s  Method  to  Coupled  Electron-Photon  Transport  jji  i he 
i  MeV  20  MeV  Energy  Range 

As  was  stated  in  the  above  section,  very  good  agreement  with  experimental  and  .Monte 
Carlo  energy'  deposition  data  was  achieved  for  electron  sources  in  the  energy  range  0.1  to  1.1 
.MeV.  in  tile  one  lugli  energy  case  that  we  treated  in  Ret.  5,  the  2U  MeV  electron  beam 
incident  on  H2O,  the  energy  deposition  profile  obtained  did  not  agree  well  with  Monte  Carlo 
results'  '  (see  Fig.  5).  It  became  apparent  that  for  high  electron  source  energies,  it  would  be 
necessary  to  include  consideration  of  radiative  energy  losses  in  our  calculations.  The  success  of 
the  coupled  electron-photon  Sn  calculations  performed  by  Lorence,  Nelson  and  Morel'  for 
isotropic  electron  sources  provided  us  with  sufficient  encouragement  to  perform  a  set  of 
calculations  for  high  energy  electron  beam  sources,  in  particular  to  repeat  the  20  MeV-IT.O 
calculation.  The  technique  developed  would  then  be  applied  to  other  beam  energies  in  H  ,0 
and  other  materials  as  well.  VVe  installed  a  bremsstrahlung  photon  production  and  transport 
module  into  our  code. 
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Energy  Deposition  (MeV*-cm2/g) 


10.0 


Beryllium,  0.1  MeV,  0  degrees 


Fraction  of  a  Mean  Range 

Figure  2.  Energy  deposition,  as  calculated  with  ONETRAN*’^  in  semi-infinite  Be  for  0.1 
MeV  electrons  incident  at  an  angle  of  0*  (Ref.  5).  Circles  (experiment)  and  histogram  (.Monte 
Carlo)  data  taken  from  Ref.  4. 


Energy  Deposition  (MeV- 


MfreiLror?rf:cid:„T«";“  “ 

(Monte  Carlo)  data  taken  from  Ref  4  ^  (experiment)  and  histogra 


Energy  Deposition  (MeV“cm2/g) 


Alumlnum/Gold/Aluminum,  1.0  MeV,  0  degrees 


0,0  0,1  0,2  0,5  0,4  0.5  0.6  0.7  0.8  0.9  1.0 

Frcction  of  o  Mean  Range 


Figure  4.  Energy  deposition,  as  calculated  with  ONETRAN  ^  ,  in  a  Al/Au/Al  configuration 
for  1.0  Me\'  electrons  incident  at  an  angle  of  0*  (Ref.  5).  Circles  (experiment)  and  histogram 
(Monte  Carlo)  data  taken  from  Ref.  4. 
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Figure  5.  Energy  deposition,  as  calculated  with  ONETRAN^'^  in  semi-infinite  HjO  for  ‘JO 
Me\'  electrons  incident  at  an  angle  of  0’  (Ref.  5).  Circles  (experiment)  and  histogram  (Monte 
Carlo)  data  taken  from  Ref.  7. 
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3.1  Elect  roil  Energ\’  Loss  C  ross  Sections 


In  our  previous  calculations  it  was  assumed  that  electron  energy  loss  was  adequately 
described  by  the  continuous  slowing  down  approximation  (CSDA).  In  ine  multigroup  Sn 
formulation,  energy  downscatter  of  electrons  was  allowed  to  occur  only  between  adjacent 
energy  groups.  This  simplifying  assumption  led  to  the  use  of  a  version  of  the  ONETRAN  code 
much  smaller  in  size  and  complexity  than  the  original.  Inclusion  of  radiative  energy  losses 
necessitated  the  use  of  the  full  O.NETRA.N  capability  to  allow  for  energy  downscatter  between 
any  two  groups.  .\s  for  the  radiative  energy  loss  cross  sections,  we  performed  our  early 
calculations  with  a  very  approximate  elementary  treatment  ,  subsequently  improved,  for  the 
bremsstrahlung  generation  cross  sections  which  proceeded  as  follows: 

the  differential  cross  section  for  the  emission  of  a  photon  with  energy  between  k(  =  ht/)  and 
k-rdk  by  an  electron  with  kinetic  energy  E  as  it  passes  through  a  scattering  medium  is  given 
by 


_  v^ToBZ^/E  +  mc,C'\ 

-JiT  -  - E - ) 


(cm"‘) 


(4) 


where  cr<,  is  a  constant, 

B  is  a  slowly  varying  function  of  E, 

Z  is  the  atomic  number  of  the  scattering  medium, 

m^c^  is  the  electron  rest  energy  (all  energies  are  expressed  in  MeV), 

N  is  the  number  of  atorns/cm^. 

The  radiative  stopping  power  (electron  radiative  energy  loss  per  unit  path  length)  is  defined  as 


=  N 


kdtr 


r'ld 


(MeV/cm) 


Substituting  Eq.(4)  into  Eq.  (5), 

1 

=  N<7oZ^(E  +  m„c^)|B  d(k/E) 
0 

=  NcroZ^(E  +  mpC*)B  , 

1 

where  B  =|b  d(k/E) 


(6l 

(7) 
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is  tlie  average  va’ue  of  B  over  the  allowed  photon  energy  range. 

If  the  parameter  B  of  Eq.(4)  is  sufficiently  slowly  varj’ing.  we  may  approximate  it  as 

B  ss  B  . 


Then  Eq.(4)  becomes 


xt<ToBZ^/'E  +  moC'\ 

-dT  ir-v"^Tr'"j  ’ 

or  4  T  ■  (cm-‘)  fOj 

Eq.  9  formed  the  basis  for  our  early  calculations  of  multigroup  bremsstrahlung  cross  sections. 
Numerical  evaluation  of  the  radiative  stopping  power  was  facilitated  for  us  by  the  fact  that  the 

(oj 

DATPAC  program  calculates  this  quantity  for  arbitrary  E  and  Z,  and  we  iiad  alrca(l\’ 
adapted  the  DATP.\C  program  to  supply  the  Mott  and  Riley  scattering  cross  sections  for 
electrons. 


The  preponderance  of  literature  pertaining  to  the  calculation  of  bremsstrahlung  cioss 
sections^ a  strong  indication  that  the  accuracy  of  our  coupled  eleclron-plnnon 
calculations  might  be  significantly  improved  by  the  use  of  a  more  sophisticated  bremsstrahlung 
cross  section  model  than  that  given  above.  By  fortunate  circumstance,  the  DAI  i'.\C 
program  incorporates  much  of  the  physics  formulated  in  the  Koch  and  Motz^**^^  paper  and 
provides  tabulations  of  the  differential  cross  section  {kd(r,.^^/dk)  as  a  function  of  eleci  ron 
energy  E  for  arbitrary  E  and  Z  values.  We  extracted  this  data  from  DATPAC.  In  Fig.  0  we 
made  a  log-log  plot  of  (d(T^i,j/dk)  vs.  k  for  the  top  ten  electron  energy  groups  in  a  100  gioup 
ONETRAN  calculation  where  we  took  the  electron  source  energy  to  be  a  20  MeV  beam 
incident  on  H2O.  Thus  the  ten  curves  correspond  to  electron  energies  ranging  from  19-9  to 
18,1  MeV  in  0.2  .MeV  increments.  The  points  indicate  the  tabulated  values  from  DATP.-XC. 
While  the  curves  are  themselves  difficult  to  read,  we  concluded  from  this  plot  and  ot!ler^ 
similar  to  it,  that  the  bremsstrahlung  cross  section  could  be  easily  approximated  by  straiglit 
lines  (on  log-log  scale)  connecting  the  points  calculated  by  DATPAC.  Thus,  our 
approximation  algorithm  for  (d<7,.aj/dk)  is  of  the  form 


~sr 


A  k*. 


(10) 


where  A  and  6  are  constant  between  any  pair  of  points  calculated  by  DATPAC. 
determination  of  A  and  b  is  made  for  each  point  pair. 
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A  unique 


(cm^/g-  MeV) 


In  oui  coupled  electron-photon  transpoit  calculation,  the  multigroup  cross  sections  for 
eiicig>  do\\  nscai ti. !  of  electron^  Li.’*  radiatue  energy  loss  were  calculated  as  follows: 
a;\ei;  ail  electio;.  ci;i'rc,>  uiuip  #;;i;c!uie  (as  si'own  in  Fig.  7i  with  G  uniloriiilv  sjiaced  electron 
energy  groujis.  eac!.  cnaractci ized  b\  midpoint  eiiergy  E,-  ^where  thei.  the  cross 

section  foi  eiiiissioi.  of  a  photon  v.ith  energy  k  i=  hr')  between  k;  and  Iw.  say.  is 


(11) 

Group 


Fig  6  Differential  bremsstrahiung  r.ross  section 
for  10  electron  energies  [19.9-18. l(  2)MeV]  as 


calculated  by  DATPAC 


Fig. 7.  Electron  energy  group  structure 


If  we  assume  that  the  energy  of  the  electron  after  photon  emission  is  E,,  then  the  limits  of  the 
integral  are 


13 


k;  =  Ev  -  E,,.. 


i  r2a) 


ko  —  Ej;  -  E,_|/2  ,  ( r2b) 

or  (E;;-  E,_[y2)’  (l'2c) 

whore  E,^;/-,  and  upper  and  lower  energy  bounds  of  group  i. 

It  may  then  be  said  that  the  dowuscalter  cross  section  from  group  g  to  group  i  due  to  tlio 
emission  of  a  plioton  witli  energy  between  k,  and  ks  is 

=  <T,.(k)  ,  (13) 

where  is  Calculated  from  Eq.s.  10  and  11  above. 

3.2  nremsstrahlunii  Photon  Source 

In  addition  to  providing  the  means  by  wdiicli  the  electron  dtiwnscatler  due  to  radiative 
energy  loss  could  be  incorporated  into  the  Sn  calculation,  the  cross  sections,  were  used 

to  determine  the  photon  source  function.  'I’he  procedure  we  followi;d  was;  1)  perform  the 
<'lcctron  transport  calculation  with  the  O .\' E'l' Il.\ .N  Sn  code;  2)  retain  a  file  of  electron  flu.\es 
a.s  functions  of  energy  (group)  and  position;  3)  calculate  the  photon  source  for  use  in  a 
s'jbseciucnt  photon  transport  calculation;  *1)  perform  the  photon  transport  calculation  with 
ONE'I  K.A.N  usiiig  mulligroup  photon  cross  sections  which  we  generated  with  the  GA.MI.EG 
code’’  5j  calculatt;  deposition  profiles  from  both  the  electron  (using  flux  file  from  step  1  and 
collision  stopping  powers)  and  photon  contributions(using  photon  flux  file  from  step  4  and 
energy  absor|rlioii  coefucients  computed  by  GAMLEG).  Tiie  energy  group  structure  u;:  the 
electron  tra.’isport  calculation  was  taken  to  be  uniform,  that  is  AE,  =  E,  -  E,^,  for  all  i.  rtio 
photon  energy  group  spacing  was  taken  to  he  uniform  for  sonic  calculation t  and  logarilhniic 
for  others.  The  transport  calculations  were  carried  out  for  lUO,  and  for  the  source  energies 
considered,  the  energy  lost  to  radiation  wa.s  generally  a  small  fraction  of  the  total  energy  loss 
by  the  electrons.  We  found  that  our  results  were  insensitive  to  the  photon  energy  group 
structure,  piissibly  due  to  the  small  overall  radiative  energy  contribution.  Also  foi  this 
reason,  v.c  dcterminc.d  that  a  second  electron  transport  calculation  arising  from  the  preseiu;e  of 
Coiniiton  and  plurto-dectrons  would  not  significjintly  alter  the  results. 
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We  made  considerable  modifications  to  the  original  version  of  the  photon  cross  section 
calculation  program  GAMLEG^*^^  The  original  version,  published  in  1965.  was  written  in 
FORTRAN  I\’  and  incorporated  a  set  of  Biggs- Lighthill  coefficients^*'"^  '^■‘compiled  in  1966.  We 
updated  the  coding  to  FORTRAN-77  and  incorporated  the  latest  available  set  of  Biggs- 
Lighthill  coefficients'  'and  appro.\imation  algorithms  for  the  cross  sections.  .A.lso  in  the 
original  version  the  photon  cross  sections  could  be  calculated  for  up  to  20  elements  per  run. 
We  changed  the  code  so  that  it  now  computes  cross  sections  for  an  arbitrarily  large  (user- 
specified)  number  of  chemical  elements.  We  configured  the  code  is  now  configured  so  that  for 
a  chemical  compound,  we  enter  as  input  the  chemical  symbol  and  weight  fraction  for  each 
constituent  element.  The  new  input  data  structure,  as  we  configured  it.  is  identical  to  that 
implemented  ir  the  ITS  code  system. 

■3. .3  Results  for  ^  MeV  and  20  Me\'  Electron  Beam  sources 

Figures  8  and  9  show  the  energy  deposition  results  obtained  with  a  100  group  coujded 
ONETRA.N  calculation  for  a  5  MeV  electron  beam  normally  incident  on  H2O.  The  "slab" 
thickness  for  the  electron  calculation  was  taken  to  be  3  cm  (the  range  of  5  MeV  electrons  in 
IIoO  is  2.G  cm.).  The  depth  of  the  scattering  medium  was  extended  to  10  cm.  for  the  photon 
calculation  to  account  for  the  backscatter  of  photons  transported,  due  to  their  long  mean  free 
path,  beyond  3  cm.  In  Figure  8  the  total  (electron  4  photo.i)  energy  deposition  profile 
calculated  with  the  ONETRAN  code  (solid  curve)  is  compared  with  a  TlGER^^^Monle  Carlo 
calculation  (histogram).  The  agreement  between  ONETRAN  and  TIGER  is  good,  with  the 
exception  of  the  region  between  1.25  and  1.8  cm.  The  electron  beam  source  energy  of  5  .MeV 
pretty  much  delineates  the  point  where  errors  due  to  numerical  straggling  become  noticable. 
We  believe  that  the  consistently  higher  ONETRAN  result  shown  between  1.25  and  1.8  cm.  is 
primarily  caused  by  numerical  straggling,  an  error  source  caused  by  using  the  multigioup 
energy  discretization  scheme  which  is  first-order  accurate.  As  will  be  seen,  this  effect  shows  up 
more  prominently  as  the  source  energy  is  increased.  In  Figure  9  we  have  plotted  the  energy 
deposition  contribution  of  the  bremsstrahlung  photons  for  the  5  MeV  electron  source  C2ise. 
The  agreement  is  generally  good,  although  the  TIGER  histogram  indicates  that  not  enough 
histories  were  run  to  achieve  statistical  reliability.  As  can  be  seen  from  the  photon  dose  scale, 
the  maximum  photon  dose  contribution  is  approximately  0.15%  of  the  maximum  electron  dose 
contribution.  For  the  electron  transport,  a  20  term  extended  transport  corrected  scattering 
cross  section  was  used  in  an  S20  ONETRAN  calculation.  For  the  photon  transport,  a 

scattering  cross  section  model  was  used  in  an  S20  ONETRAN  calculatio.n. 
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Figure  8.  Energy  deposition  (electron  +  photon)  profile,  as  calculated  with  ONETRAN^*' 

f  2 1 

(curve)  and  with  ITS-TIGER  (histogram),  in  semi-infinite  HjO  for  5  MeV  electrons  incident 
at  an  angle  of  0*. 


IG 
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Figure  9.  Energy  depoeition  (photon)  profile,  as  calculated  with  ONETRaN^*^ 

(curve)  and  with  ITS-TIGER*’^(hUtogram),  in  semi-infinite  HjO  for  5  MeV  electrons  incident 
at  an  angle  of  0*. 
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In  Figures  10  and  11  are  plotted  the  total  (electron  +  photon)  and  photon  energy 
deposition  profiles  obtained  with  a  dO  group  coupled  ONETRAN  calculation  for  a  20  MeV 
electron  beam  noiuially  incident  on  H;0,  12  cm.  thickness  (20  Me\’  electron  range  in  H-^O  is 
9.3  cm.).  .-\s  in  the  previously  described  calculation,  an  extended  scattering  medium  depth  was 
assumed  (_GU  cm.;.  In  Figure  10  the  energy  deposition  profile  calculated  with  tlie  O.NETH.A.N 
code  (solid  curve)  is  compared  with  a  TIGER  Monte  Carlo  calculation  (solid  histogram)  and 
an  EGS4^‘^  Monte  Carlo  calculation  (dashed  histogram).  The  agreement  between  ONETRAN 
and  TIGER  is  excellent  to  a  depth  of  about  5.5  cm.,  beyond  which  the  O.N'ETRAN  result  is 
consistently  higher.  The  EGS4  calculation  includes  knock-on  electrons  which  account  for  the 
lower  energy  deposition  in  the  first  1  to  1.5  cm.  (when  knock-ons  are  included  in  the  TIGER 
calculation,  the  two  .Monte  Carlo  results  agree  very  well).  The  photon  dose  results  (Fig.  11) 
agree  fairly  well.  W’c  believe  that  the  discrepancies  found  may  be  caused  in  part  by  our  crude 
model  for  the  photon  source  angular  distribution.  The  effects  of  numerical  straggling  on  the 
electron  energy  deposition  profile  (Fig.  10)  are  definitely  more  pronounced  for  the  20  Me\’ 
source  than  for  the  5  MeV  case.  This  was  somewhat  mitigated  by  increasing  the  number  of 
energy  groups,  but  even  when  we  increased  the  number  of  groups  to  100,  the  straggling  effect 
persisted  and  was  significant. 

MK  nock-on  Elect ron  Production 

In  modelling  high  energy  electron  transport,  we  made  modifications  to  our  electron  cross 
section  module  that  would  include  electron  energy  losses  and  secondary  electron  production 
through  knock-on  collisions.  As  was  stated,  the  EGS4  calculation  shown  in  Fig. 10  included  the 
effects  of  knock-on  electron  production,  while  TIGER  was  run  for  the  case  of  no  knock-ons. 
W'lieii  knock-ons  were  taken  into  account,  the  dose  near  the  surface  of  incidence  was  focud  tu 
be  about  -5  percent  lower  than  that  obtained  without  the  knock-ons.  Evidently,  the  highest 
energy  electrons  suffc.ed  a  sufficient  amount  of  knock-on  collisions  to  experience  significant 
downscatter  in  energy.  As  primary  electrons  penetrate  the  medium,  the  continuous  slowing- 
down  collisional  and  radiative  energy  losses  outweigh  any  visible  effect  of  downscatter  rhrough 
knock-on  collisions.  In  order  to  resolve  this  disagreement  between  Monte  Carlo  and  discrete 
ordinates,  we  instalUed  a  knock-on  algorithm  in  the  calculation  of  the  electron  energy  loss 
cross  sections  for  the  ONETRAN  calculation. 
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Figure  10.  Energy  deposition  (electron  4  photon)  profile,  w  ceiculnted  with  ONETRAN 
(curve)  with  ITS-TlGER'*^(*olid  histogram)  and  EGS4^^\dasbed  histogram),  in  semi-infinite 
HjO  for  20  MeV  electrons  incident  at  an  angle  of  O’. 
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Figu  re  11.  Energy  deposition  (photon)  profile,  as  calculated  with  ONETRAN^'^ 

f2l 

(curve)  and  with  ITS-TIGER  ^histogram),  in  semi-infinite  H2O  for  20  MeV  electrons  incident 
at  an  angle  of  0*. 
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Accounting  for  downscatter  of  primary  electrons  and  secondary  electron  production  due  to 
knock-on  collisions,  translated  into  the  calculation  of  multigroup  cross  sections  for  these 
proceses.  Along  with  this,  it  was  found  necessary  to  use  a  modified  stopping  power  function  to 
account  for  the  fact  that  some  of  the  energy  of  the  primary  electrons,  rather  than  being 
entirely  e.xpended  along  the  primarj-  track,  i.s  transfered  to  secondary  electrons,  wliich  are  in 
turn  downscaitered.  The  definition  of  “restricted  collision  stopping  power"  and  the  means  for 
its  calculation  are  provided  in  Ref.  16. 

Our  calculation  of  the  multigroup  knock-on  cross  sections  proceeded  as  follows: 
the  total  knock-on  cross  section  (trj'’)  for  primary  electrons  in  group  g  with  energy  E,  is  given 
by 

^  tUda  i  =  1  /  -  J 

=  f  'i-') 


where 


(15) 


and 


(E,u,  0.01  MeV) 


(16) 


E,<e  and  a:  •  respectively,  the  secondary  electron  energy  and  low  energy  cutoff  for 

the  secondary  electrons, a  ad^*^' 


d<T  _  C  /  1  ,  1  ,  /  r  V  (2r-fl)  /  _1_\1 

^  -  I3^E,\  P  (1-c)*  \r+lj  (7+17  ^ 


(!') 


with  C  =  2)rNar<,mc^^^^  (  =  .153536^^^),  r  =  (  =  the  group  energy  in  inc- 

units).  /?  is  the  electron  velocity  divided  by  the  speed  of  light,  T.  o  =  Avogadro’s  number,  r9  = 
electron  radius,  Z  and  A  are  the  atomic  number  and  weight,  respectively,  and  signifies 
the  average  Z/A  for  the  molecular  composition  of  the  scattering  material. 

If  we  allow  the  maximum  fractional  energy  loss  of  the  primary  to  be  Cmax  =1/2.  as 
indicated  in  the  upper  limit  of  the  integral  (Eq.l4),  then  the  total  knock-on  cross  section  for 
group  g  is 
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.V  nt,  -  Krrfei)}- 


In  view  of  the  above,  and  given  the  electron  energy  group  structure,  uniformly  spaced 
energy  groups  of  width  AE,  as  depicted  in  Fig,  12,  the  cross  .section  for  and  electron  with 
energy  Ej,  to  produce  a  knock-on  with  energy  E^,  and  therefore  suffer  a  reduction  in  its  own 
energy  to  E^„  (=  E^  -  E^,  )  is 


^ 

.g 


ii^y 


(  ^ g' /  ^  ^marj 


(  ^ g' y  —  ^mar  ) 


where 


t,  =  mm 


fj  =  max 


e^,+.1ae 

Eg 

i 


s  tmat/i, 


Combining  Eqs.  17  ,19-21  we  obtained  for  the  group-to-group  cross  section, 


O'  fit 


^  /  .  —if  _  ^  _ r:y _  ^  _ . 

/Jy-Ej  E  ,+Ue  E  ,-iAE  E»-E  ,.iAE  E»-E  ,  +  1aE 


+  .  (iim  i„f  ( 


-y  f  +  tr  V  '  E  ,  -  iAE''  E,-E  ,-iAE  ^  ) ) 


Eq.  22  does  not  include  the  special  cases  where  <i  or  t-g  —  (.max.  These  endpoint  cases  wen- 
treated  in  our  computer  program. 

After  a  knock-on  collision  has  occurred,  the  primary  electron,  with  energy  E^„  has  undergone 
angular  deflection  v  given  by^'^^ 


Energy 


Group 


f 


The  angle  x  between  the  knock-on 
secondary  and  the  incident  primary 
electron  directions  is  given 

The  Legendre  coefficients  of  the  angular 
deflection  dependence  of  the  primary 
electron  dovvnscattei  cross  section 
a  ,,  are  then  given  by 

g  —  g 

It 

(t/”'  =  ^,,P;(cOs-.)  . 

The  corresponding  Legendre  coefficients 
for  the  direction  of  the  secondary 
electron  are 

<T^_^,^P,(cosx)  .  (23) 


With  some  experimentation,  we  found  that  (max  could  be  adjusted  to  impose  realistic 
limits  on  the  primary  electron  energy  loss  due  to  knock-on  collisions.  While  the  theoretical 
value  for  (nax  =  1/2,  stepwise  examination  of  the  integrated  values  of  Eq.l7  revealed  that 
such  occurrences  were  rare.  The  best  agreement  with  Monte  Carlo  was  obtained  with  the 
following  upper  limits  on  the  knock-on  electron  energy:  for  primary  energy  5  MeV  and  above, 
2  MeV;  for  primary  energies  below  5MeV,  200  keV.  To  account  for  the 

energy  carried  off  by  the  secondary  electrons,  a  separate  transport  calculation  for  the  knock- 
ons  was  performed.  Figure  13  shows  the  total  energy  deposition  profile,  obtained  using  the 
knock-on  algorithm,  for  the  20  MeV  beam  source  on  HjO.  Even  though  agreement  with  the 
EGS4  result  is  much  improved  in  the  first  1.5  cm.  over  that  shown  in  Fig.  10,  this  calculation, 
also  in  contrast  to  the  Fig.  10  curve,  consistently  underestimates  the  dose  in  the  region  beyond 
the  profile  maximum.  It  is  clear,  from  the  above  results,  that  improvement  in  both  the 
knock-on  algorithm  and  the  energy  discretization  s,  >  .le  is  required. 


Fig  12  Electron  energy  group  structure  showing  relative 
positions  of  primary  and  secondary  electron  energies 
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Figure  13.  Energy  deposition  (electron  +  photon)  profile,  as  calculated  with  ONETRAN^  ^ 
(curve)  with  ITS-TIGER'^'(solid  histogram)  and  EGS4^^^(dashed  histogram),  in  semi-infinite 
HjO  for  20  MeV  electrons  incident  at  an  angle  of  O’.  The  knock-on  electron  contribution  is 
included  in  this  calculation. 


M  Structure  of  ih£  Coupled  Electron-Photon  Discrete  Ordinates  Programs 


We  wrote  and/or  modified  programs  for  the  one-dimensional  coupled  electron-photon 
discrete  ordinates  calculation  with  a  modular  design  consisting  of  four  sections: 

1)  an  interactive  (with  the  user)  program  to  collect  the  user-specified  information  which 
defines  the  problem  to  be  run  and  translates  this  information  into  an  input  data  file 
formatted  to  be  compatible  with  the  cross  section  and  transport  codes. 

2)  cross  section  generation  programs  which  are  modified  versions  of  DATP.\C^'’  (for 
electrons)  and  GAMLEG^*^\for  photons). 

3)  data  handling  programs  to  couple  the  electron  and  photon  transport  calculations  by 
converting  and  passing  data  between  them;  compute  quantities  such  as  energy 
deposition  profiles,  energy  and  angle  spectra;  prepare  plot  data  files. 

4)  O.N’ETRAN^*^  transport  code. 

Of  the  above,  the  programs  in  the  second  and  fourth  categories  are  well  documented  by 
their  originators.  The  linking  programs,  third  category,  perform  perfunctory  programming 
tasks.  The  first  program,  the  user-calculation  interface,  which  is  run  interactively,  has  self- 
explanatory  prompts,  and  the  input  format  from  the  keyboard  is  free  form.  We  wrote  two 
versions  of  this  program,  one  in  BASIC  and  one  in  FORTRAN-77,  so  that  it  can  be  run  either 
on  an  IBM  PC  or  compatible  (BASIC  version)  or  on  a  mainframe  such  as  VAX  or  MicroV.^X. 
The  user  is  interrogated  regarding  such  problem-defining  items  as  the  number  of  mat  ''r.ls  in 
the  scattering  medium;  the  chemical  composition  (element  symbol  and  weight  f.ac*>.  .i  of  <»ach 
material;  source  type  (electron,  photon,  angular  distribution,  spatial  distrubu'.  a.  eiw  / 
spectrum);  whether  a  photon  transport  calculation  is  to  be  run;  whether  a  knock-cn  c/.culation 
is  to  be  included;  etc.  This  program  also  creates  the  list  of  files  (with  definitions  and  names 
linked  in  a  directory)  to  serve  as  repositories  for  data  (such  as  electron  and  photon  fluxes, 
cross  sections,  plot  data).  These  files  are  used  by  the  data  handling  programs  (category  3, 
above)  and  ONETRAN.  This  program  also  contains  provisions  for  eliminating  redundant 
calculations  of  cross  section  sets  by  providing  the  user  with  the  opportunity  to  specify  the 
names  of  already  existing  cross  section  files  which  may  be  pertinent  to  the  problem  to  be  run. 
A  run  of  this  program  results  in  the  creation  of  an  “input  run  deck”  for  the  transport 
calculation  on  the  mainframe.  If  the  BASIC  version  is  used,  the  “run  deck”  file  is  written  to  a 
floppy  diskette  which  can  then  be  transferred  to  the  mainframe  for  a  batch  mode  transport 
calculation.  A  hard  copy  of  the  prompts,  responses,  run  file  and  a  data  file  directory  ic 
printed. 
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ilL  ELECTRgN  TRANSPORT  CALCULATIONS  IN  ONE-.  TWO-  AND  THREE 
DIMENSIONS 
Introduction 

During  this  reporting  period,  there  were  many  opportunities  to  provide  benchmark 
electron  transport  calculations  for  the  purposes  of  testing  newly  developed  one-  and  two- 
dimension  electron  transport  techniques  and  scattering  cross  section  formulations.  Most  of  the 
benchmark  calculations  were  made  using  the  Monte  Carlo  method.  In  addition  to  performing 

fl9] 

benchmark  Monte  Carlo  and  Sn  calculations  in  one-dimension  ,  we  also  wrote  a  new  single 

[21] 

collision,  electron  Monte  Carlo  (1-d)  transport  code  for  the  purpose  of  testing  the 
effectiveness  of  a  new  scattering  cross  section,  the  “smart”  cross  section. 

Our  participation  in  research  in  the  area  of  two-dimensional  electron  transport  calculation 

methods  consisted  of  providing  benchmark  Monte  Carlo  calculations,  using  the  ITS  Serie-s 

f22] 

ACCEPT  code,  for  two  papers.  The  first  '  demonstrated  the  feasibility  of  using  the  "sniart" 

cross  section  in  conjunction  with  first  collision  (or  once-scattered)  sources  for  electron  beam  Sn 

[231 

calculations.  The  second  of  these  *  'represents  a  continuation  of  the  electron  beam  source 
work.  It  contains,  however,  a  mere  extensive  set  of  comparisons  between  Sn  and  Monte  C  arlo 
results  consisting  of  isodose  contour  maps  and  energy  deposition  profiles  in  one-dimension 
across  the  two-dimensional  scattering  medium. 

We  conclude  this  chapter  with  a  discussion  of  ihrce-dimeiis'onal  electro.’i  tran-sj-ort 
calculations.  A  recent  paper'*^'  which  reports  the  results  of  a  set  of  Sn  and  Monte  Carlo 
calculations  of  energy  deposition  in  and  electron  exit  currents  from  realistic  microelectronic 
device  geometries  due  to  the  incidence  of  a  beam  of  electrons.  We  provided  the  device 
geometry  description  and  mock-up  as  we'l  as  the  Monte  Carlo  calculations  for  this  paper. 

Monte  Carlo  and  Sn  Calculations  of  Electron  Transport  in  Otie-Dimension  for  Validation  ^ 
the  Metliod  of  St  ream  inn  Ravs 

In  a  paper  which  appeared  in  1987^*'*^  a  complete  description  was  given  of  the  application 
of  the  method  of  streaming  rays  to  one-dimensional  electron  transport  problems.  Earlier 
papers  {op.  cit.)  had  already  demonstrated  the  applicability  of  this  calculation  method  to 
neutral  particle  transport  in  two-dimensional  media  with  void  regions  and  time-dependent  one¬ 
dimensional  problems.  The  similarity  bet  .veen  the  time-dependent  neutral  particle  transjujit 

26 


equation  and  the  Spencer-Lewis  equation  for  electron  transport  strongly  suggested  the 
applicability  of  the  streaming  ray  (SR)  method  to  this  problem.  In  the  Spencer-Lewis 
equation,  the  electron  energy  is  the  continuous-slowing-down-approximation  (CSDA)  defines  a 
unique  relationship  between  the  electron  pathlength  and  stopping  power.  The  reason  for 
employing  this  approximation  is  that  the  major  mechanism  for  electron  energy  loss,  inelastic 
collisions  with  other  electrons,  cannot  be  treated  in  a  tractable  way  by  a  form  of  the  transport 
equation  that  considers  one  collision  at  a  time.  There  are  far  too  many  of  these  collisions,  and 
the  electron  energy  loss  resulting  from  an  individual  inelastic  collision  is  miniscule.  The  SR 
method  is  highly  compatible  with  CSDA  because  the  SR  model  assumes  that  the  motion  of  an 
electron  is  assumed  to  take  place  along  a  streaming  ray  (a  characteristic  line  of  the  Spencer- 
Lewis  equation).  With  the  SR  method  the  electron  pathlength/energy  increases/decreases 
continuously.  This  is  the  basic  point  of  departure  from  the  conventional  multigroup  Sn 
method  in  which  the  spatial  dependence  of  the  Spencer- Lewis  equation  solution  is  obtained  by 
finite  differencing  over  a  series  of  discrete  steps  for  a  given  constant  energy  value.  When,  by 
means  of  an  iterative  spatial  differencing  scheme,  the  equation  is  balanced  for  a  particular 
energy  group,  the  electron  energy  is  then  decremented  by  a  discrete  step  to  the  next  lo<.ver(in 
energy)  group. 

Numerical  errors  associated  with  the  spatial  and  energy  discretization  of  the  Sr>  metliod 
give  rise  to  artificial  straggling.  In  some  cases  this  leads  to  a  fortuitously  close  agreement  with 
Monte  Carlo  and  experimental  measurements  of  energy  deposition  because  the  Monte  Carlo 
method  is  capable  of  accouinting  for  the  actual  physical  straggling,  the  departure  from  CSD.\. 
Numerical  straggling  in  the  Sn  method  introduces  large  errors  in  the  calculated  electron  energy 
spectra  unless  a  very  large  number  of  energy  groups  is  used.  This  'effect  is  greatly  reduced 
with  the  SR  method.  Ref.  19  provides  a  method  for  simulating  physical  straggling  in  addition 
to  CSDA  in  the  transport  calculation.  Fig.  14  shows  a  comparison  of  energy  deposition  results 
obtained  with  the  SR,  and  Monte  Carlo^^^  methods.  These  results  were  obtained  for  the 
case  of  a  slant  source  (po  =  .916667  ,  see  Fig.  14  inset)  of  0.2  MeV  electrons  on  a  0.04  g/cm^ 
A1  slab.  The  SR  calculation  was  formulated  and  performed  by  Filippone,  while  the  Sn  and 
Monte  Carlo  calculations  were  done  by  Woolf.  As  can  be  seen,  the  SR  calculation,  which  has 
straggling  effects  included,  agrees  better  with  the  Monte  Carlo  benchmark  than  either  of  the  2 
Sn  calculations. 
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Figure  1-4.  Energy  depuiition  profile  for  0.2  MeV  electrons  obliquely  incident  (pg  =  0.916C67) 
oil  ft  0.0-1  g/cni'  thick  Al  slab  (Ref.  19). 
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Figure  15.  CSD.A  transmission  spectra,  as  calculated  by  the  SR,  Monte  Carlo  and  So  methods, 
for  0.2  MeV  electrons  obliquely  incident  (po  =  0.916667)  on  a  0.04  g/cm®  thick  Al  slab  (Ref, 
19). 
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Fig.  15  sliows  the  CSD.A  transmission  spectra  for  tlie  same  0.2  Me\'  source  geometry. 
However  in  both  the  SR  and  Monte  Carlo  calculations.  CSD.\  was  assumed  i,that  is,  the 
straggling  effects  were  omitted  on  purpose).  The  "couittershove"  correction,  referred  to  in  the 
figure,  is  a  method  devised  by  Filippone  to  correct  the  spatial  discretization  errors  associated 
with  cell-collided  electrons.  SR  treats  the  collided  and  uncollided  portions  of  the  electron  flux 
separately.  The  uncollided  part  is  treated  exactly,  while  the  collided  part,  generally  much 
smaller  if.  for  example,  extended  transport  corrected  cross  sections  are  used,  is  treated  in  the 
same  manner  as  in  Sn.  As  can  be  seen  from  Fig.  15.  the  SR  and  Monte  Carlo  results  aggree 
very  well,  while  even  the  SO  group  Su  calculation  is  far  off  the  mark  above  0.13  MeV.  The 
delta  function  referred  to  in  the  figure  occurs  at  the  maximum  energy  at  which  the  uncollided 
electrons  can  emerge  fiom  the  slab,  given  the  slant  path  and  CSDA.  Fig.  16  compares  the  SR 
and  Monte  Carlo  benchmark  calculations  with  the  inclusion  of  physical  straggling.  Again  the 
agreement  between  the  two  methods  is  good. 


E  <MeV) 


Figure  16.  Transmission  spectra,  as  calculated  by  the  SR  and  Monte  Carlo  methods,  for  0.2 
MeV  electrons  obliquely  incident  (/tq  =  0.916667)  on  a  0.04  g/cm*  thick  A1  slab.  Straggling 
effects  are  included  in  both  calculations  (Ref.  19). 
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^  Imnlementation  of  SMART^^°*  Scattering  Theory  in  One-Diinensional  Monte  Carlo  Electron 

(2i] 

Transport  Calculations 

f  '>f  1 

The  SMART'*  scattering  cross  section,  originated  by  Filippone,  is  a  representation  of  an 
angular  redistribution  function  for  electron  scattering  based  on  Goudstnit-Saunderson  theory 
(see  Ref.  20).  The  Goudsmit-Saunderson  multiple  scattering  theory  allows  one  to  obtain 
infinite  medium  solutions  for  the  Spencer-Lewis  equation.  The  basis  of  SMART  scattering 
theory  is  the  conjecture  that  the  validity  of  the  scattering  theory,  and  its  representation 
through  the  scattering  matrix  should  not  depend  on  the  problem  geometry.  A  rigorous 
treatment  of  SMART  scattering  theory  is  given  in  Rei.  20,  thus  it  is  not  necessary  or 
appropriate  to  repeat  it  here.  .\s  a  result  of  collaboration  between  Drs.  Filppone  (Univ.  of 
Arizona)  and  Woolf  (Arcon  Corp.)  a  paper  was  written  and  published,  first  in  a  meeting 
proceedings'^*"^,  and  then  in  Nuclear  Science  and  Enymeeriny*^*^^describing  the 
implementation  of  this  cross  section,  in  scattering  matrix  form,  in  a  Monte  Carlo  electron 
transport  code  written  by  Dr.  Woolf.  The  term  "SMART”  is  a  acronym  for  simulating  many 
accumulative  Rutherford  trajectories.  These  scattering  matrices  were  originally  developed  for 
implementation  with  the  discrete  ordinates  (Sn)  method.  One  of  the  purposes  of  performing 
these  Monte  carlo  calculations  was  to  demonstrate  that  the  validity  of  the  SMART  scattering 
theory  did  not  depend  on  the  method  of  implementation;  that  its  applicability  to  the  Sa 
method  wa«  not  due  to  a  numerical  artifact  of  the  method.  The  other  purpose  of  this  work 

was  to  introduce  a  possible  alternative  means  for  performing  electron  Monte  Carlo  calculations 

[2] 

to  the  implementation  of  multiple  scattering  theory  as  manifested  in  the  ITS  code  series  .  It 
was  felt  that  application  of  SMART  scattering  theory  in  the  Monte  carlo  method  could  result 
in  a  more  economical  means  for  performing  these  calculations  than  those  which  presently  exist. 
This  is  due  in  large  part  to  the  fact  that  an  essential  characteristic  of  this  theory  is  a  large 
effective  mean  free  path  for  electrons,  much  large  than  the  true  single  collision  mean  free  path. 
The  extent  to  which  this  feature  would  result  in  a  more  economical  calculation  than  that 
afforded  by,  say,  the  TIGER^^^code,  which  employs  a  multiple  collision  cross  section  and  is  also 
much  more  economical  than  true  single  collision  Monte  Carlo,  remained  to  be  determined  by 
the  results  of  this  paper. 

The  Monte  Carlo  code  written  for  the  implementation  of  the  SMART  scattering  theory, 

incorporated  a  single  collision  trajectory  analog  algorithm.  It  departed  from  conventional 

Monte  Carlo  algorithms  in  two  important  ways;  1)  the  calculation  of  the  intercollision  distance 
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takes  into  account  the  variation  of  the  total  cross  section  with  electron  path  length  and 
direction:  and  2)  the  direction  cosines  of  the  electron  trajectory  after  a  collision  has  occurred 
are  obtained  directly  from  the  scattering  matrix.  The  total  electron  trajectory  path  length, 
determined  from  the  electron  CSDA  range,  was  divided  into  a  number  of  equal  path  length 
steps  As.  The  intercollision  path  length  was  determined  by  .‘am|)ling  the  well-known 
exponential  attenuation  laws.  When,  along  the  intercoliision  parth  segment,  a  given  As  step 
was  traversed,  the  total  cross  section,  a  function  of  path  length  and  direction,  was  updated. 

The  presently  existing  form  of  the  Monte  Carlo  code  treated  one-dimensional  slab 
geometry  problems.  Extension  to  two-  and  three-dimensions  would  be  straightforward.  For 
purposes  of  testing  the  code,  two  200  keV  electron  source  configurations,  beam  and  isotropic, 
were  assumed  to  be  incident  on  two  thicknesses  of  Al  ,  .01  g/gm'  and  .02g/cm'.  The  incident 
obliquity  of  the  beam  source  was  defined  by  its  cosine  Hq—  0.916667.  The  results  obtained 
with  this  code,  emergent  energy  spectra  and  energy  deposition  profiles,  were  compared  with 
benchmark  runs  of  TIGER  .  A  sampling  of  the  results  given  in  the  paper  are  shown  in  Figs 
17-20.  A  comparison  of  the  transmission  spectra  (SMART  and  TIGER)  through  0.01  g/cm' 
Al  is  shown  in  Figs.  17  and  18  for  the  beam  and  isotropic  sources,  respectively.  Beam  source 
reflection  spectra  for  the  thicker  slab  are  compared  in  Fig.  19.  Finally,  the  energy  deposition 
profiles  for  the  beam  source  in  .02  g/cm*  are  compared  in  Fig.  20.  As  can  be  seen,  the  results 
obtained  with  the  SMART  Monte  Carlo  calculation  compare  very  well  with  the  TIGER 
benchmarks.  Both  codes  were  run  with  10000  electron  case  histories.  In  all  cases  tested,  it 
was  found  that  the  SMART  matrix  Monte  Carlo  calculations  ran  approximately  3  times  faster 
than  TIGER.  Since  our  code  was  developed  for  research  rather  than  production  purposes,  it  is 
likely  that  this  run  time  ratio  could  be  improved. 
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Figure  19.  Reflection  spectrum  for  a  0.2  Me\'  electron  beam  source,  with  incident  obliquity 
=  0.916667,  incident  on  a  0.02  g/cni*  .\1  slab,  as  calculated  by  SM.\RT  cross-section  Monte 
Carlo  and  ITS-TIGER'^'  .Monte  Carlo  (Ref.  21). 
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Figure  20.  Energy  deposition  profile  for  a  0.2  MeV  electron  beam  source,  with  incident 
obliquity  /tq  =  0.916667,  incident  on  a  0.02  g/cm’  A1  slab,  calculated  by  SMART  cross- 
section  Monte  Carlo  and  ITS-TIGER  ^  Monte  Carlo  (Ref.  21). 
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4.  Electron  Transport  Calculations  in  Two- Dimensions 

4.1  Determination  of  Dose  Profiles  and  Isodose  Contours  for  200  keV  Electron  Transport  in 
Two-Dimensions 


In  a  paper  that  resulted  from  a  collaboration  among  W.  L.  Filippone,  S.  Woolf  and  J.C. 
Garth  *  ,  it  was  reported  that  the  Spencer-Lewis  equation  for  electron  transport  in  two- 
dimensions  was  solved  using  the  Sn  method  with  diamond  differencing  in  the  spatial  and 
pathlength  variables.  Dose  profile  calculations  in  two-dimensions  were  made  for  200  keV 
electrons  incident  on  Al.  The  electron  source  configuration  were  of  three  types;  1)  normally 
incident  beam;  2)  beam  with  slant  incidence  (45*);  and  3)  isotropic  incidence.  Comparison 
runs  were  also  made  with  the  ACCEPT  ‘  Monte  Carlo  code. 


In  two-dimensions,  the  Spencer-Lewis  equation  for  the  electron  flux  <^(x,y,s,n)  is  given 


{§s  +'^(s))  <?(x,y,s,n)  =  I  dnV(s,fi'-*n)(^(x,y,s,n')  +  Q(x,y,s,fi)  ,  (24) 

4t 

where  Dx,  Hy  =  the  x  and  y  velocity  direction  cosines, 
s  =  electron  pathlength, 
cr(s)  =  total  scattering  cross  section. 
cr(s,il^-*f2)  =  differential  scatering  cross  section, 
and  Q(x,y,s,f])  =  electron  source  density. 

The  energy  deposition  function,  W(x,y),  is  determined  from  the  Spencer-Lewis  solution 
4>{\,y,s,Q)  as  follows: 

W(x,y)  =  ||^l<?i(x,y,s)  ds,  (25) 

where 


1^1  is  the  electron  stopping  power,  and  <^(x,y,s)  is  the  scalar  flux. 

Eq.  24  is  solved  using  the  Sn  method  with  diamond  differencing.  The  solution  steps  are 

outlined  in  Ref.  22.  Because  of  the  extreme  anisotropy  of  the  electron  scattering  kernel, 

SMART^^'^'  scattering  theory  was  employed  to  render  the  Sn  numerical  solution  feasible.  An 

effective  scattering  matrix  was  defined  in  order  that  the  required  number  of  discrete  ordinates 
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for  the  description  of  the  angular  redistribution  of  electrons  would  not  exceed  manageable 
limits.  Figures  21  and  22  are  isodose  contour  plots  of  the  energy  deposition  in  Al  for  normal 
and  45‘  electron  beams,  respectively.  In  both  cases,  the  dimensions  of  the  Al  rectangle  were 
0.01  X  0.02  g/cm'.  For  a  direct  quantitative  comparison  with  the  Monte  Carlo  method,  we 

r  T  1 

ran  the  ,\CCEPT  code,  which  is  a  three-dimensional  electron  Monte  Carlo  program,  for  the 
two  cases  shown  in  Figs.  21  and  22.  The  two-dimensional  problem  geometry  was  simulated 
by  extending  tne  z-coordinate  of  the  scattering  medium  to  essentially  ±oc,  where  “cc”  is  taken 
to  be  ~1.5  electron  range  units.  Figures  23  and  24  show  comparisons  of  dose  profiles  obtained 
with  Sn  calculations  and  ACCEPT.  These  profiles  are  the  integrated  (over  area)  closes 
across  the  rows  A.  B.  and  C  as  indicated  in  the  accompanying  cell  diagram.  In  Fig.  23.  the 
isotropic  source  results  for  two  Sn  calculations,  Sg  and  Si2»  are  compared  with  Monte  Carlo, 
while  in  Fig.  24,  the  corresponding  comparisons  for  the  normal  beam  source  results  for  S.,  and 
Sg  is  are  shown.  .‘\s  can  be  seen,  overall  agreement  between  our  Monte  Carlo  calculations  and 
Filippone’s  S,,  calculations  is  very  good. 
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4.2  T’wQ-DifTiensional  Electron  Transport  Calculations  with  Qnce-Collided  Beam  Sources 


In  a  continuing  collaboration  among  W.L.  Filippone,  S.  Woplf  ap/d  J.C.  Garth  on  two- 

(23] 

dimensional  electron  transport  problems,  a  paper  was  presented  at  the  November.  1987 
meeting  of  the  American  Nuclear  Society  in  which  the  use  of  analytic  first  collision  (or  once- 
collided)  electron  sources  was  introduced  into  two-dimensional  Sn  electron  transport 
calculations.  In  previous  one-dimensional  electron  transport  calculations,  it  was  found  that  Sn 
calculations  could  not  be  done  for  beam  sources  unless  an  analytically  derived  once-scattered 
source  was  used.  This  was  due  to  the  high  degree  of  anisotropy  associated  with  electron 
scattering  kernels.  While  in  two-dimensional  electron  transport  Sn  calculations  incorporating 
the  SMART  cross  section,  this  anisotropy  could  be  handled  without  the  aid  of  the  once- 

Mol 

scattered  source  ,  it  was  believed  that  use  of  such  a  device,  ,iu  conjunction  with  the  SMART 
cross  section,  could  greatly  improve  the  accuracy  of  these  calculations.  As  in  the  paper 
discussed  in  section  4.1,  we  provided  the  Monte  Carlo  benchmark  calculation  for  this  work.  A 
200  keV  beam  of  electrons  was  assumed  incident  on  an  A1  slab  of  dimension  .01  g/cm'  x  .02 
g/cm^.  The  angle  of  incidence  in  the  x-y  plane  was  45*,  and  the  point  of  incidence  was  located 
at  x=0,  y=.01,  as  shown  on  the  inset  of  Fig,  25.  Also  shown  in  Fig.  25  are  the  two  energy 
deposition  results,  Sn  (  S4)  and  Monte  Carlo  (ACCEPT^*',  lOOOOO  case  histories).  To  aid  in 
the  interpretation  of  the  diagrams  of  Fig.  25,  it  should  be  pointed  out  here  that  the 
rectangular  slab  was  divided  into  50  square  zones  (.002  x  .002),  and  the  histograms  display  the 
energy  (keV)  deposited  in  each  zone.  The  results  appear  to  agree  well.  The  advantage  in 
computation  speed  of  Sn  over  Monte  Carlo  became  especially  apparent  in  this  work.  It  was 
found  that  the  S^  calculation  consumed  ~l/30-th  as  much  computer  time  as  did  the  Monte 
Cailo. 
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2i  FJectron  Transport  iri  Three-Dimensions 


Throughout  all  of  the  previous  one-  and  two-dimensional  electron  transport  research  in 
which  we  have  participated,  a  primary  objective,  in  addition  to  achieving  high  computational 
accuracy,  has  been  to  achieve  a  high  degree  of  realism  in  our  modelling  capability.  In  pursuit 
of  the  first  objective,  consideration  of  realistic  three-dimensional  device  geometries  was 
postponed  until  the  present.  Also  causing  delay  in  this  regard  was  the  unavailability,  until 
recently,  of  a  good  three-dimensional  electron  Monte  Carlo  program,  and  finally,  the  lack  of 
adequate  de.scriptions  of  realistic  microelectronic  device  geometries.  The  availability  of  the 
ITS^^'  Monte  Carle  code  series  has  satisfied  the  first  requirement,  while  the  appearance  of  two 
papers*^*  "^'  in  the  1987  IEEE  Transactions  on  Nuclear  Science,  in  which  two  microelectronic 
device  geometries  were  published,  satisfied  the  second.  The  work  to  be  discussed  here,  resulted 
in  the  presentation  of  a  paper  at  the  American  Nuclear  Society  Mathematics  and  Computation 
Topical  Meeting^^^'.  This  paper  will  appear  in  Nuclear  Science  and  Engineering  in  early  1990. 

In  Ref.  24  is  reported  the  extension  of  Filippone’s  Sn  method  for  solving  the  Spencer-Lewis 
equation  for  electron  transport  to  treat  three-dimensional  multiregion  problems.  One  of  the 
key  points  of  this  work  was  the  utilization  of  the  flux  continuity  condition,  which  is  generalized 
for  muitiregion  problems  by  expressing  the  flux  as  a  function  of  electron  energy  rather  than 
pathlength.  Such  a  representation  results  in  a  set  of  Sn/diamond  difference  equations  which 
are  nearly  identical  in  form  to  conventional  Sn/diamond  difference  equations.  The  Sn  method 
was  then  used  to  calculate  electron  energy  deposition  due  to  200  keV  electron  beams  incident 
on  problem  geometries  typical  of  silicon  and  gallium-arsenide  semiconductor  microelectronic 
devices.  Our  roles  in  this  work  were:  1)  to  research  and  define  the  realistic  device  geometries; 
2)  to  provide  quantitative  descriptions  of  these  silicon  and  gallium-arsenide  device  geometries 
for  modelling  purposes;  3)  to  define  the  inputs  required  by  the  ACCEPT  Monte  Carlo 
calculation;  and  4)  to  aid  in  the  interpretation  and  presentation  of  the  results. 

A  fairly  extensive  discussion  of  the  theory  of  the  three-dimensional  Sn  calculation  method 

was  included  by  Filippone  in  Ref.  24,  and  thus  need  not  be  reviewed  here.  The  result  of  the 

electron  transport  calculations  were  presented  in  the  form  of  energy  deposition  resulting  from 

the  incidence  of  electron  beam  sources  on  problem  geometries  typical  of  semiconductor 

microelectronic  devices.  The  configurations  we  chose  to  work  with  are  based  on  fairly  detailed 

descriptions  of:  1)  a  silicon  (Si)  hybrid  diode-resistor  and  2)  a  gallium-arsenide(GaAs)  field- 

effect  transistor^^^^.  The  final  configurations  for  the  device  geometries  that  were  used  in  our 
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calculations  are  simplified  versions  of  these,  however  their  essential  geometrical  character  and 
material  compositions  were  preserved.  The  Si  device  geometry  is  shown  in  Figure  26. 

The  Sn  code  was  run  for  200  keV'  electron  beam  sources  incident  at  a  slightly  off-center 
point:  x=0yrm,  y=2S/im,  z=:0/«m  for  the  Si  device  (as  shown  in  Fig.  25).  The  incident  beam 
orientation  is  given  by  ^=90®,  ^=45®,  where  6  is  the  polar  obliquity  with  respect  to  the  z-axis 
and  <i>  is  the  azimuthal  angle  in  the  x-y  plane.  The  scattering  cross  sections  were  taken  to  be 
screened- Rutherford. 

VVe  made  independent  calculations  of  energy  deposition  using  the  ACCEPT  module  of  the 
Integrated  TIGER  jeries(ITS)^*W  coupled  electron-photon  Monte  Carlo  codes.  The  ACCEPT 
code  was  run  using  the  screened- Rutherford  cross  section  with  the  CSDA  (no-straggling) 
option.  In  this  way  it  was  possible  to  ensure  that  both  the  Sj^  and  Monte  Carlo  calculations 
were  based  on  the  same  physical  model.  Sufficient  numbers  of  Monte  Carlo  case  histories 
(10®)  were  run  so  that  valid  benchmark  results  could  be  obtained.  The  estimated  standard 
error  was  <1  percent  in  all  of  the  Monte  Carlo  tally  cells  for  which  the  energy  deposition 
exceeded  1  percent  of  the  maximum  value.  The  worst  cases,  5-6  percent  error  estimates,  were 
encountered  in  only  2  of  129  cells  for  the  Si  calculation  and  1  of  120  cells  for  the  GaAs 
calculation.  In  these  instances  the  energy  deposition  values  were  less  than  0.1  percent  of  the 
peak  value. 

The  energy  deposition  results,  Sn  and  Monte  Carlo,  for  the  Si  device  mock-up  are  shown  in 
Figures  27a  and  27b,  respectively.  The  data  shown  result  from  integration,  over  the  z- 
coordinate,  of  the  energy  deposited  per  unit  volume  in  the  region  bounded  by  the  gold  cap,  and 
are  therefore  given  in  units  of  MeV/cm*.  The  energj'  deposition  results  obtained  with  Sn  were 
found  to  compare  well  with  results  of  ACCEPT  Monte  Carlo  calculations.  The  agreement 
between  the  Sn  and  Monte  Carlo  results  for  the  GaAs  device  configuration  is  qualitatively 
similar.  This  comparison  is  shown  in  Ref.  24.  Computer  run  times  required  for  the  Sn 
calculations  were  found  to  be  lower  than  that  required  for  Monte  Carlo  by  factors  ranging 
from  30  to  50. 
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Figure  26.  Silicon  device  geometry  (not  drawn  to  scale).  Lenglli  units  are  /im  (Ref.  2d) 


Figure  27a.Three-dimensional  Sn  calculation  of  energy  deposition  (MeV/cm^)  resulting  from  a 
•15*  electron  beam  incident  in  a  Si  device.  Source  energy  =  0.2  MeV  (Ref.  24). 


Figure  27b.ITS-ACCEPT*^^  Monte  Carlo  calculation  of  energy  deposition  (MeV /cm*)  resulting 
from  a  45*  electron  beam  incident  in  a  Si  device.  Source  energy  =  0.2  MeV  (Ref.  24). 


J V,  LOW  ENERGY  ELECTRON  TRANSPORT  (fl  -  2fl  fiY) 


1.  Introduction 


We  have  written  computer  programs  and  performed  calculations  to  model  the  inelastic 
scatter  and  transport  of  electrons  with  very  low  energy  (0  -  20  eV)  in  SiOj.  The  inelastic 
scatter  can  occur  via  two  processes:  1)  ionization  and/or  excitation  collisions  with  other 
electrons;  ^.nd  2)  electron-phonon  interactions.  The  first  of  these  can  occur  in  SiOo  only  for 
electron  energies  above  the  ionization  threshold  of  8.9  eV.  It  is  believed  that  the  second  of 
these  processes,  interactions  with  phonons,  can  occur  across  the  energy  range  0  -  20  eV.  Our 
principal  activities  in  this  research  area  have  consisted  of  three  projects.  The  first  of  these  is 
tlie  interpretation  of  electron-phonon  scattering  theory  and  its  translation  into  a  set  of 
numerical  algorithms,  and  a  computer  program  (PHONON),  for  the  evaluation  of  basic 
scattering  parameters,  such  as  scattering  cross  section,  scattering  angle  distributions,  and 
effective  stopping  powers.  The  second  task  in  this  project  is  the  incorporation  of  the  above 
calculated  parameters  in  a  Monte  Carlo  transport  code  ^LOWEND)  for  electrons  in  the  U  -  20 
eV  energy  range  in  SiOj-  The  third  of  these  is  the  application  of  the  method  of  discrete 
ordinates  (Sr,),  a  determininstic  method  for  the  numerical  solution  of  the  transport  equation, 
to  electron-phonon  scattering. 

Zi  PHONON:  Computer  Program  for  the  Evaluation  af  lil£  Basic  Parameters  of  Electron- 
Phonon  Scattering 

Electrons  of  very  low  energy,  below  ionization  threshold,  can  undergo  inelastic  scatter 

through  the  electron-phonon  interaction.  These  are  interactions  with  the  lattice  and  result  m 

phonon  emission  or  absorption.  Specific  information  relating  to  the  evaluation  of  electron- 

phonon  cross  sections,  stopping  powers  and  scattering  angles  comes  to  us  by  way  of  the  both 

private  communication  and  open  literature  .  The  theoretical  basis  of  the  contents  of  this 

section  is  adapted  from  the  work  of  Ashley  and  applies  to  electron-phonon  scattering  in 

SiOj-  With  this  theoretical  btisis,  we  designed  the  numerical  algorithms  and  wrote  the 

computer  program  PHONON  for  the  evaluation  of  the  cross  sections,  stopping  powers  and 

scattering  angle.  Some  of  the  results  of  our  calculations  to  be  shown  here  also  appear  in  Ret. 

28.  Electron-phonon  interactions  occur  in  two  modes,  longitudinal  (LO)  and  acoustic  (ACj 

There  is  a  total  of  6  interaction  modes,  4  LO  and  2  AC.  This  arises  from  the  fact  that  both 

the  LO  and  AC  modes  may  result  in  phonon  emission  (electron  energy  gain)  or  phonon 
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absorption  (electron  energy  loss),  and  there  are  2  discrete  phonon  energies  to  consider  for  Si02 
(hwp  =:  0.068  eV  and  0.153  eV[wp  is  the  phonon  frequency]). 

2.1  LO  Modes 


2.1.1  Inverse  Mean  Free  Patli 


The  inverse  mean  free  paths  for  the  LO  modes  are  given  by  the  following  expression^ 


where  <^7  LO  ~  inverse  mean  free  path  for  phonon  emission, 

A2  LO  ~  inverse  mean  free  path  for  phonon  absorption, 
wp  =  phonon  frequency, 

v  =  incident  electron  velocity, 

<00  =  Si02  dielectric  constant  in  the  optical  region, 

<0  =:  Si02  static  dielectric  constant, 

nw  =  ^exp(/lw/kBT)  -l)"‘. 

kgT  =  .025  eV, 

q  =  wave  number  corresponding  to  electron  energy  gain  (+)  or  loss  (-)  . 


The  total  inverse  mean  free  path  (cross  section)  for  the  LO  mode  is  then  given  by 

-^LO  “  '^+,LO  +  (2^) 

The  upper  and  lower  limits  of  the  integral  in  Eq.  25  correspond,  respectively,  to  the  maximum 
and  minimum  possible  q  values:  if  the  incident  electron  kinetic  energy  is  denoted  as  E,  then 

Wj±(E)  =  min{ko(l+^lT^  )  .  ^bz  }  .  (28a) 

and 
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W±(E)  =  ±ko(l-JlT^  ), 


(28b) 


where 


and 


ftko  =  ^2m*E  , 

ni'  =  electron  effective  mass, 

^BZ  wavenumber  corresponding  to  the  Brillouin  zone  threshold  energy. 


The  quantity  -^]  has  two  values,  0.0816  for  ftwp=0.068  eV  and  0.112  for  flu;p=0.153  e\'. 


2.1.2  Scattering  Angle  Distribution 


The  distribution  of  scattering  angle  cosines,  cos^^,  is  determined  from  momentum  and 


enerf^y  conservation  and  is  given  by 


[27] 


cose  (E)  =  l-qV2kg±ftu»p/2E 

cos^^  ,q(L)  _ 


<291 


where,  as  before,  the  momentum  tranfer  wave  number,  q,  ranges  from  Wj^(E)  to  VV^(E). 


The  average  scattering  angle  cosine  <  cos0^  ^_q(E)>  is  obtained  from 


W±(E) 


<  cos^^  =  I  ^(l-qV2kg±ftwp/2E)  / flT^Wp/E 


Wj±(E) 


q 


W±(E) 


W±(E) 


(lTftwp/2E)  ln[W±(E)/W±(E)]  -  ^2(Wi*=(E)'-W±(E)'] 
(,jlq:ftwp/E  )  ln[V/2=^(E)/W*(E)r  ~ 


VoU; 


2.1.3  Stopping  Power 


The  expression  for  the  LO  mode  stopping  power  is  given  by 


[27] 


dE| 

dx 


(nw+i±^)ln 


l.n.-Wj*=(E) 

W±(E) 


(31) 
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2.2  AC  Modes 


2.2. 1  Inverse  Mean  F ree  Path 

The  inverse  mean  free  paths  for  the  AC  modes  are  given  by 


qmaz 


132) 


where  M  and  N£  are,  respectively,  the  mass  and  density  of  the  SiOo  unit  cell 
and 

T  -111  C$/fi  ,  q  <  ^B2 


r  'Jko  ?  2n\  C,, 

I  ko  +  JlTftC,k 


+  A  lTftCfk02/E^  ,  q  >kQ2 


(33) 


{  (e^  Vi)-1 


,  q  <  kQ2 
’  ^  -^B2 


(34) 


r  C.q 

u/p(q) 

I  Cik02 


*  <l<ltBZ 


C<  =  velocity  of  sound  in  Si02,  and 


a  =  ftC,/k0T 


(3o) 


2.2.2  Scattering  Angle  Distribution 


For  the  AC  modes  the  distribution  of  scattering  angle  cosines  cosd^  is  also  determined 

l57] 

from  momentum  and  energy  conservation  and  is  given  by 

*■*'='  jl=Fllu,(q)/E 

where  the  phonon  frequency  u'p(q)  is  now  dependent  on  the  the  momentum  transfer  wave 

number,  q,  as  given  in  Eq.  35.  In  contrast  to  the  LO  case,  an  exact  analytical  expression 

cannot  be  derived  for  the  average  scattering  angle  cosine  <  cos^^  ,ac(E)>- 
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„  ± 
^max 


<  cos(?^  ^^(E)>  = 


^  [n(q)  +  i±ll(l-qV2kg±ft--,/2E) 
dq — - - 


0 


j  qmax 

/I 


dd  ^  [n(q)+.^±.^(3r) 


must  be  evaluated  numerically. 


2.2  -3  Stopping  Power 

The  expression  for  the  AC  mode  stopping  power  is  given  by 


[27] 


qmax 


dE| 

dlT'AC^ 


=  3m' 

4jrMNc/lkfl 


I  dq  q®  [n(q)  + 


(38) 


2.3  N'umcrical  Evaluation  of  Scattering  Parameters 

The  above  expressions  for  the  electron-phonon  inverse  means  free  paths,  scattering  angle 
distributions  and  stopping  powers  were  evaluated  in  our  computer  program,  PHOXO.N,  Since 
the  expressions  for  the  longitudinal  mode  parameters  (given  by  Eqs  26,30,31)  are  analytic  and 
do  not  involve  complicated  numerical  evaluation  of  special  functions,  tables  of  these  quantities 
as  functions  of  electron  energy  were  generated  in  a  straightforward  manner.  Algorithm 
development  for  numerical  evaluation  of  the  acoustic  mode  parameters  was  somewhat  more 
involved.  The  integral  expressions  (Eqs.  32,37,38)  could  not  be  evaluated  analytically  witliout 
approximation.  We  divided  the  region  of  integration  into  3  sub-regions:  1)  0  to  q,,  where  the 
energy  fraction  ftw^/kgT  =  0.005  :  2)  q;  to  kQ2  ;  and  3)  kg^  qmax  . 

In  the  fir.st  sub-region,  recalling  Eq.  35,  we  made  the  following  approximations  : 


(  1  :f  ftC.q/E)~‘^^  «  1  ±hC,q/2E  , 


and 


.  V  kgT 


(39) 

(40) 


With  these  approximations,  the  integrands  of  Eqs.  32,37,38  become  simple  algebraic  functions 
of  q  which  are  then  evaluated  exactly. 


The  second  integration  sub-region,  q,  to  kg^,  is  handled  with  our  adaptive  Gaussian 
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quadrature  numerical  integration  program,  AQG4.  This  program  employs  the  method  of 
interval  splitting  to  perform  the  numerical  integration  to  within  a  user  specified  accuracy.  If 
such  accuracy  is  unattainable,  due  to  the  functional  form  of  the  integrand  or  computer 
limitation.  .\QG4  informs  the  calling  program  through  informative  diagnostics. 

Finally,  in  the  third  sub-region,  both  --pCq)  and  n(q)  take  on  constant  values,  so  that  the 
integrals  of  Eqs.  32,37,38  are  evaluated  analytically.  We  made  extensive  calculations,  for  the  4 
LO  and  2  AC  interaction  modes,  based  on  the  above  theory  for  the  electron-phonon  scattering 
cross  section,  mean  scattering  angle  cosine  and  stopping  power  in  SiOi  for  subionization 
electrons  (E  <  8.9  eV).  In  Fig.  28  we  show  a  plot  of  the  total  scattering  cross  section, 
averaged  over  the  LO  and  AC  modes.  Fig.  29  is  a  plot  of  the  mean  scattering  angle  cosine 
averaged  over  the  LO  and  .^C  phonon  emission  modes  (solid  curve).  The  dashed  curves  show 
the  contributions  on  the  LO  and  AC  emission  components  individually.  Finally,  Fig.  30  is  a 
plot  of  the  electron-phonon  interaction  stopping  power,  averaged  over  the  6  modes. 

Tabulations  of  the  electron-phonon  cross  section,  average  scattering  angle  cosine  and 
stopping  power  were  compiled  using  the  PHONON  program  and  then  incorporated  into  ojr 
low  energy  electron  transport  Monte  Carlo  trajectory  simulation  program,  LOWEND  (to  be 
discussed  in  the  next  section).  It  was  then  possible  to  obtain,  through  trajectory  simulation, 
electron  energy  spectra,  energy  loss(gain),  and  drift  velocity  values  in  the  presence  of  an 
applied  electric  field  for  subionization  electrons  in  Si02- 


ScolUrIng  Crest  Section  (Angstroms-) 


0.40  - 


Figure  28.  Total  electron  scattering  cross 
section,  a\eraged  over  the  longitudinal 
optical  and  acou.'tical  phonon  excitation 
inodes,  in  SiO^.  for  the  sub-ionizatiou 
threshold  energy  region  (Kef.  28). 


Figure  29.  Mean  scattering  angle  cosine  for 
electrons  with  energy  <  8.89  eV  in  SiO;  for 
longitudinal  optical  phonon  (LOE)  and 
acoustical  phonon  (ACE)  emission.  The  solid 
curve  is  the  weighted  mean  cosine  over  the 
two  phonon  excitation  modes  (Ref.  28). 


Figure  30.  Electron  stopping  power 
averaged  over  the  longitudinal  optical  and 
acoustical  pl.oncn  excitation  modes,  in 
SiOj  for  the  sub-ionization  threshold 
energy  region  (Ref.  28). 
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^  LQVVEND:  Mg^t^  Carlo  Simulation  Program  Subipnizat[9.n  (<8.9  £\Q  Electron 
Transport 


3.1  Overview 

The  physical  parameters  governing  the  transport  process  for  low  energy  (subionization) 
electrons  were  calculated  with  program  PHONON  and  then  supplied  through  a  set  of  tables  to 
our  single  scatter  Monte  Carlo  program  LOWEND'  These  tables  consisted  of  the  following 
energy-dependent  entries  for  the  electron-phonon  interaction:  1)  total  inverse  mean  free  path, 
A"‘;  2)  total  stopping  power  ,  dE/dx:  3)scattering  angle  average  cosine,  <cos-/>.  The  total 
inverse  mean  free  path  and  stopping  power  values  were  obtained  by  summation  over  tlie  4 
longitudinal  and  2  acoustic  phonon  emission  and  absorption  modes,  while  the  scattering  angle 
cosine  values  were  obtained  with  a  summation  over  the  6  modes  weighted  by  their  respective 
probabilities  of  occurrence. 


The  main  features  of  the  LOWEND  code  are  that  it  accounts  for:  1)  electron  energy  loss 

and  scatter  with  both  LO  and  AC  phonons;  and  2)  energy  loss/gain  in  the  presence  of  an 

applied  electric  field.  With  this  code  it  is  possible  to  obtain,  through  trajectory  simulation. 

electron  energy  spectra,  energy  loss(gain)  or  A£  distributions,  and  drift  velocity  values  in  the 

presence  of  an  applied  electric  field  for  subionization  electrons  which  appear  in  SiOj  as  a  result 

of  inelastic  scattering  processes.  The  electric  field  values  for  which  this  program  has  been 

tested  range  from  0  to  10  MV/cm.  In  this  calculation  electron  trajectories  are  simulated  in 

direct  analog  Monte  Carlo  fashion.  It  was  found  that  5000  -  10000  electron  case  histories  were 

sufficient  to  produce  statistically  adequate  results  for  energy  spectra,  while  1000  histories 

proved  sufficient  for  determination  of  averaged  quantities  such  as  drift  velocity.  The  above 

estimates  appied  provided  that  the  case  histories  were  allowed  to  run  for  a  sufficiently  large 

number  of  simulated  collisions  to  produce  quasi-equilibrium  conditions  (small  variation  of 

o 

energy  spectrum  shape  and  magnitude)  over  a  large  spatial  region  (300-1000  A).  This 
condition  gave  rise  to  a  “rule  of  thumb”  that  600-1000  collisions  per  case  history  were 
adequate  for  electric  field  (S)  values  below  5  MV/cm,  while  a  larger  number,  -^3000  for 
determination  of  averaged  quantities  and  '>'10000  for  energy  spectra,  were  required  for  6  =  10 
MV/cm.  Typical  run  times  for  the  Monte  Carlo  program  aie  on  the  order  of  0.6  hr/(1000 
histories)  on  the  AFGL  VAX  8650  computer. 
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3.2  Electron  Trajectories  in  the  Presence  of  an  Applied  Electric  Field 

The  action  of  the  applied  field.  8.  is  to  transform  the  electron  trajectory  path  segments 
into  parabolic  sections.  For  short  path  segments,  this  curvature  can  be  ignored  and  the 
calculations  can  be  made  on  the  basis  of  a  straight  line  segment  approximation.  The  length  of 
the  segment,  As,  is  chosen  by  random  selection  from  the  path  length  probability  distribution 
(via  the  usual  inversion  of  the  exponential  attenuation  rule,  A(E)  is  taken  to  be  the  mean  free 
path).  The  energj  change  is  then  given  by 

As/Vg 

AE  =  I  F  ds  =  I  F-ydt 


As/Vg 

=  £-(vc 


V  +  at)dt  . 


%» here 


HP  HP 

Fj  =  ^  V  +  eSz,  ^  =  electron  stopping  power. 


a  =  V  +  ||z.  e  =  electron  charge. 


V'Z  =  cos6  =  velocity  direction  cosine  along  8  di.ection. 


Insertion  and  integration  yields 

AE  =  (f  +  .ecos())As  +  (|JE|>  +  (e6)>  +  2  .  («) 

From  these  expressions  one  can  calculate  the  velocity  components  and  position  coordumtes. 

For  longer  As  segments,  which  admit  of  curvature  and  larger  energy  change,  we  liavi- 
used  a  five  sub-segment  approximation  to  the  path  integral  for  energy  change.  In  this 
procedure  the  values  of  E,  cos6,  A(E),  dE/dx  and  the  position  and  velocity  coordinates  are 
updated  in  each  step.  The  resulting  energy  change  is  given  by 


Ej  -  F,  =  +  eScos^i)  + 

k  k  -ux  / 


-f  (e8)2  +  2  ^e8cos5*)(^~  -  . 
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The  velocity  component  (v^^)  in  the  direction  of  the  electric  field  (8)  is  given  by 


of  Electron  Energy  Spectra 


One  of  the  results  obtainable  with  a  LOWEND  calculation  is  the  determination  of  kinetic 

energy  equilibrium  spectra  in  the  presence  of  an  applied  electric  field.  When  we  considc'red 

electron  transport  exclusively  in  the  subionization  energy  region,  it  was  found  that  for  a  •*<  e\' 

o 

injection  energy  into  SiOo.  energy  equilibrium  was  well  established  at  a  depth  of  500  A  for  a 
wide  range  of  electric  field  values.  Fig.  31  shows  calculated  electron  kinetic  energy  spectra  for 
3  values  of  8  ,  (0.5,  4  and  10)xl0®V/cm.  When  we  extended  the  electron-phonon  cross 
section,  scattering  angle  and  stopping  power  tables  to  the  8.9  -  20  eV  range, ionizing  collisions 
occurred.  The  ionization  cross  sections'  '  used  in  this  calculation,  supplied  by  ORNL,  are 
shown  in  Fig.  32.  Fig.  33  shows  the  electron  energy  spectrum  in  the  presence  of  a  lO'V/cm 
electric  field.  This  result  was  obtained  using  electron-phonon  parameter  tables  extended  to  20 
eV  and  the  OR.N'L  inelastic  cross  sections.  As  can  be  seen,  the  high  energy  portion  of  the 
spectruin(>S.9  eV)  accounts  for  a  significant  fraction.  It  was  found  that  for  fields  in  excess  of 
'>'7xlO®V/cm.,  our  predictions  of  the  amount  of  electron  multiplication  exceeded  experimental 
observation.  An  effort  is  presently  underway  to  reconcile  these  differences. 
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Figure  33.  Kinetic  energy  spectrum,  as  calculated  by  LOWEND,  for  applied  field  of  lO^V/cm 
at  500  A  depth  in  SiOj.  High  energy  tail  above  8.9  eV  exists  due  to  small  cross  section  for 
secondary  ionization  at  these  energies. 
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4.  Transport  Equation  Solution  for  Electron-Phonon  Scatterini^  bv  the  Method  of  Discrete 
Ordinates  (Sn) 

Discrete  ordinates  electron  transport  calculations  were  made  using  the  standard  S„  code, 
ONETRAN^^',  for  low  energy  (subionization)  electrons  in  SiOj.  Preparation  for  these 
calculations  involved  reformulating  the  electron-phonon  transport  parameters,  cross-section 
and  stopping  power  into  the  multigroup  format  compatible  with  the  3n  method.  As  was  done 
for  the  Monte  Carlo  calculations  reported  in  the  previous  section,  these  multigroup  parameters 
were  derived  separately  for  the  longitudinal(LO)  and  acoustic(AC)  mode. 

4.1  LQ  Modes 

4.1.1  Group-to-Groun  Scatterimr  CiToss  Sections 


Making  use  of  Eqs.  26,  28a, b  and  29,  we  arrived  at  expressions  for  the  group-to-gioup 
energy  transfer  cross  sections. 

For  energy  loss: 

let  Eg  =  initial  electron  energy,  and  E^'  =  final  electron  energy  ,  so  that  Eg>Eg'  . 
Then  the  downscatter  cross  section  is 


-E/)  =  I 


Aw  n,.  ln[q^^,/Wr(E,)]  ;  Wr(E,)  <  <  Wj-(E,) 

I 

0  ;  elsewhere 


where 


=  ,j2m’(E,-E,')  . 


There  are  two  cross  sections  —'Eg')  corresponding  to  the  two  values  of 

Similarly,  for  E^<Ej'  (energy  upscatter)  we  have 


rtO(E,  -*E/)  =1 


Aw  (nu,+l)  ln(q^^,/Wr(E,)]  ;  Wr(E,)  <  q^^,  <  W^(E,) 
0  ;  elsewhere 


where  ^^gg'  ~  >l2m’(Ej'-Ej) 

As  in  the  i-honon  absorption  case,  there  are  two  cross  sections  -»Ey')  corresponding 

to  the  two  values  of  ftw. 
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4.1.2  Group-to-Group  Scattering  Angle 


The  scattering  angle  cosines  corresponding  to  the  above  are  given  by 

1-q]  ./2kUftW2E, 

- JitL/e.  • 

where  tikg  =  ^2m*Ey  . 

Expressions  were  derived  for  the  Legendre  expansion  coefficients  of  the  angular  scattering 
cross  section  for  electron-phonon  interactions.  This  was  necessitated  by  the  fact  that  most 
discrete  ordinates  algorithms  and  codes  accept  and  utilize  the  angular  scattering  data  in  this 
format.  The  analysis  for  the  longitudinal  optical  (LO)  mode  downscatter  in  energy  from  group 
g  to  group  g'  {Ey>Ey')  is  shown  below  (the  derivation  for  upscatter  was  similarly  derived). 


The  average  scattering  angle  cosine  is  given  in  terms  of  integrals  over  the  momentum 
transfer,  q,  by; 


w* 


where 


I  ^  (l-qV2k5+flw/2E,) 


</i-  />  = 

39 


w-  , 

1.99 


'^0  ' 
i.99 


jTWE;  j 

W" 


dq 

■?r 


1.99 


'^7  I  =  k  j2m*(E,-E,Vi/j)  , 

l.gg  n 


W-  ,  =  i  .|^■(E,-Ey'-l/2)  , 

7,93 


with  the  restrictions  that 


and 


minIw-^^,]=  Wr(E,)=  k,(l-|I+^  )  , 
max[w-^  ,]  =  WjCEy)  =  min|k,^l-f  ^  ),  kB2}  • 


(48) 


(49a) 

(49b) 


(50a) 

(50b) 


where  hk 


^2m'E  ,  m*  is  the  electron  effective  mass  and  kg^  tke  wavenumber 


corresponding  to  the  Brillouin  zone  threshold  energy. 
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Tlie  Legendre  expjansior.  for  the  angular  scattering  cross  section  is 


=  ]^(/-fl/2)cr,P;(A<)  . 


♦  i 

<T,  =  [fftpjP,(^)  d/r  .  (52)  Energy 

^1/7  - 

El - 

With  the  following  definition  for  the  n-th  moment  of  ,  E3/2  “ 

E.. - 


Group 


o{^^)  p''  d// 


(r{^t)  dp 


the  first  five  Legendre  expansion  coefficients  are  given  as 


+  i 

= 


(54a) 


Ej-1/2  - 

E, - 


<Tj  =  <n>ffQ  ,  (54b) 

<f-2  =  ^(3<p->  -  l)<ro,  (54c) 

(Ti  =  J(5</i^>  -  3<p>)ro,  (54d) 

=  i(35<p^>  -  30<p’>  +3)<r(,.  (54e) 


Figure  34.  Energy  group  structure 
for  subionization  electrons,  SiOj 
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The  above  moments  of  ^  were  calculated  v/iih  analytical  integrations  of  the  form 


I  ^  (l-qV2k|+ftw/2E,)'‘ 

- w*'  ,  — —  •  (55) 

S9  2,ig' 

il-'i-'/Ey  I  ^ 

Program  PHONON  (see  sec.  2  of  this  chapter),  which  was  originally  written  to  calculate 
electron-phonon  cross  sections  for  the  LOWEND  (see  .sec.  3  of  this  chapter)  Monte  Carlo 
calculations,  was  modified  to  incorporate  the  above  algorithms  for  the  Legendre  expansion 
coefficients. 

AC  Modes 

The  acoustic  mode  (AC)  cross  sections  were  treated  differently  from  the  LO  cross 
sections.  The  simplifying  assumption  of  isotropic  scattering  was  made  for  the  angular 
dependence  of  the  AC  cross  sections.  This  arises  from  the  following  considerations.  The 
maximum  energy  transfer  in  an  AC  collision  is  ~0.02  eV.  The  average  scattering  angle  cosine 
for  these  individual  collisions  is  ~-0,3  and  is  fairly  constant  until  the  energy  approaches  values 
less  than  2  eV.  For  a  reasonably  sized  energy  group  width,  say  0.1  eV,  outscatter  to  other 
groups  would  not  occur  in  any  single  collision.  It  is  known  that  many  such  individual  collisions 
can  occur  (for  much  of  the  energy  range  of  interest  Thus  it  was  decided  to 

treat  the  AC  mode  case  using  the  continuous  slowing-down  approximation.  Since  we  already 
had  expressions  for  the  AC  mode  stopping  power  (up-  and  down-8catter)as  given  by  Eq.  38, 
we  cou'd  calculate  the  cross  section  for  transfer  between  adjacent  energy  groups  as 

‘^AC,,  =  ’  (56) 

where  **  fhe  AC  mode  stopping  power. 

The  above  cross  section  (Eq.  56)  may  be  considered  as  representing  a  composite  of  many 
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individual  collisions,  the  number  of  which  depends  on  the  energy  group  width  AE.  This 
composite  ‘•collision"  consists  of  many  turn-around  collisions  which  average  out  to  a  scattering 
angle  cosine  0  (isotropic  scatter).  This  algorithm  was  also  incorporated  in  program 
PHONON.  In  addition,  a  section  of  code  was  wriiten  for  PHONON  to  convert  the  cross 
section  Legendre  coefficient  table  into  a  format  compatible  for  use  with  ONETRAN. 

4.3  Results  of  Sn  Calculations 

A  series  of  preliminary  electron  transport  calculations  were  made  with  O.NETRAN  for 
electron-phonon  scattering.  The  calculations  were  done  in  the  absence  of  an  electric  field.  The 
electron  source  energy  was  taken  to  be  5  eV,  The  calculations  were  made  with  100  uniformly 
spaced  energy  groups  over  the  energy  range  0  to  10  eV.  The  thicl  ^ss  of  the  scattering 
material,  SiOj,  was  taken  to  be  0.1  pm,  with  100  spatial  discretization  intervals  of  equal  width 
(Az  =  O.OOi  pm).  Four  source  geometries  were  considered:  1)  planar(at  z  =  Oj 
monodirectional  (along  z-direction)  source  [fig.35a];  2)  monodirectional  (along  z)  source 
uniformly  distributed  from  z  =  0.04  pm  to  z  =  0.06  pm  [fig.  35b];  3)  monodirectional  (along  z) 
source  of  width  Az  centered  at  z  =  0.0505  pm  [fig.  35c];  4)  isotropic  (in  direction)  source 
uniformly  distributed  from  z  =  0.04  pm  to  z  =  0.06  pm  [fig.  35d]. 


Z-* 

0  O.lpm 

Fig.  .35a. 

0  O.lpm 

Fig.  35b. 

Fig.  35c. 


0.1pm 


Fig.  35d. 


O.lprn 


Sample  plots  of  thz  electron  flux  profiles  as  calculated  for  5  energies  are  shown  in  Figs.  36 
and  37  for  the  plane  perpendicular  and  isotropic  source  configurations  shown  in  Figs.  35a  and 
35d,  respectively. 
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Plane  Monodirectionol  Source  ot  Z 
Source  Energy  =  5  eV 


9|oi:^jDd  sojnos/xn|j  uojp9|3 
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Figure  36.  Electron  (lux  vs.  depth  prordes.  as  calriilatvil  with  ONF/l'IlAN  *  for  electron- 
phonon  scattering,  for  5  electron  energies,  due  to  the  iiicidriice  of  .1  5  oV  nioiiodircction.tl 
electron  source  normally  iiichicnt  on  0.1  SiOj. 


4.4  Treatment  of  Electron  Transport  in  tiig  Presence  sf  aa  External  Electric  Field 


Calculations  of  electron  transport  in  the  presence  of  an  applied  external  electric  field  have 
been  done  by  Wienke^"^^  using  the  method  of  discrete  ordinates.  We  began  the  process  of 
modifying  our  O.N'ETRAN  calculations  to  include  an  electric  field  term.  We  used  as  a 
theoretical  basis  the  method  outlined  in  Ref.  29. 

The  Lorentz  force  term,  ff^(E  +'v'xB)  Vv^,  is  incorpcrated  into  the  right-hand-side 
(source  side)  of  the  Boltzmann  equation.  The  quantities  e,  me,  E  and  Bare,  respectively, 
the  election  charge,  mass  velocity,  electric  and  magnetic  fields.  ^  is  the  electron  flux.  For  our 
purposes,  the  magnetic  field  portion  was  omitted.  In  this  formulation,  electric  fields  are 
treated  as  effective  sources.  The  electric  fields  are  regarded  as  having  the  effect  of  up-  and 
down-scattering  electrons  in  energy  and  also  redistributing  their  directions.  Slowing  down  and 
speeding  up  of  electrons  due  to  electric  fields  are  treated  in  the  continuous 
approximation(CSD A),  i.e.,  only  energy  transfers  between  adjacent  groups  are  allowed. 
Assuming  slab  geometry,  Wienke’s  treatment  of  the  Lorentz  force  source  term  (electric  field 
portion  only)  proceeded  as  follows: 

defining  the  direction  cosine  p  in  one-dimensional  plane  geometry, 

T-T  IS  vx/i,  (57) 

th"  gradient  operator  (velocity  space)  is  given  as 

^•'=-^^  +  (1^+^)^,  (58) 

one  obtains  (n^)E‘^v^  =  (  fT^)bE^-l-^(l-p*)^]  .  (59) 

The  first  term  in  [  ]  above  is  the  energy  upscatter  or  downscatter  term,  depending  on  the 
direction  of  E,  while  the  second  term  in  [  ]  is  the  angular  redistribution  term.  The  first  term 
can  be  expressed  in  finite  difference  form  as 

(  ’  for  /i  >  0  (60) 


and 
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(  for  /i  <  0  ,  (Gi) 

where,  as  in  Fig.  34,  Ey  >  Ej+,  and  <  =  (  fj^E). 

The  subscript  m  denotes  tiie  direction  of  the  m-th  discrete  ordinate,  nm.  Both  cases  above 
contribute  a  positive  downscattered  or  upscattered  source  to  group  g  in  the  direction  m.  The 
upscatter  terms  necessitate  outer  (energy)  iterations  in  the  Sn  calculation. 

The  angular  redistribution  term,  was  treated  as  a  within-group  collision  (self-scatter) 
term.  Since  in  Sn  calculations,  the  angular  dependence  of  the  flux  is  expressed  in  terms  of  a 
Legendre  series  expansion  of  order  M-1  in  p,  it  is  advantageous  to  make  use  of  the  fact  that 

in  the  flux  expansion,  so  that 


(  -£L)  (1 V)^  -  £  (2/+l/4tr)f' , 


,0  _  e  r(2M^-2M-l),  ,  _  o  ,  r(2/^  •  21-1), 

’  ~  ''^H2M  +  l)(2M-3)^  '  ’7^1(2M)(2/-|-3)J  ‘ 


(2M)(2/-t-3)^ 


The  extended  transport  correction  of  order  M-1  was  used  in  the  electric  field  expansion. 

When  the  above  expressions  were  incorporated  into  the  discrete  ordinates  form  of  the 
transport  equation,  all  of  the  terms  containing  4>$m  are  then  transposed  to  the  left-hand  side 
(removal  term8)of  the  transport  equation  and  grouped  into  the  total  loss  term  ,  with 

the  appropriate  modification  of  the  total  cross  section  (<ry).  The  remaining  source  terms 
involving  group-to-group  transfer  and  angular  redistribution  will  be  grouped  on  the  right-hand 
side  (gain  terms)of  the  transport  equation  with  the  in-group  source,  the  fixed  source, 
secondary  electron  source  and  inelastic  downscatter  source. 

Since  Wienke  states  in  Ref.  29  that  the  analysis  described  above  can  be  used  to  modify 

existing  Sn  codes,  we  began  the  process  of  incorporation  of  the  electric  field  terms  into  our  Sn 

calculations.  This  work  is  presently  in  progress.  We  believe  that  the  use  of  another  Sn  code, 
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rather  than  ONETRAN,  may  be  preferable.  We  have  a  code  which,  unlike  ONETRAN, 
directly  calculates  electron  fluxes  (rather  than  flux  Legendre  moments),  and  makes  direct  use 
of  the  source  terms  (rather  than  their  Legendre  representation).  Such  a  code  may  prove  more 
suitable  for  the  incorporation  of  electric  field  terms  than  ONETRAN. 
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